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COMPUTER AND ENGINE PERFORMANCE STUDY OF A GENERALIZED 
PARAMETER FUEL CONTROL FOR JET ENGINES 
by Kurt Seldner and Harold Gold 
Lewis Research Center 

SUMMARY 

A mathematical analysis of a generalized parameter hydraulic fuel- control concept 
is presented. An analog computer simulation was used to establish the feasibility of the 
fuel-control concept for turbojet engine applications. The simulation of the fuel control 
was first operated with a simulation of the J85-13 engine and then operated as an experi- 
mental control with an actual J85-13 engine in a test cell. Results obtained from the use 
of the simulated fuel control with both the simulated and actual engines are presented. 
The operation of the control is discussed, and its performance is compared with that of 
the normal J85-13 control. 


INTRODUCTION 

An analysis of a fuel-control concept for turbojet engines is presented in this report. 
The analysis results from a study conducted at Lewis as a portion of a program to inves- 
tigate engineering problems that affect production costs of turbojet engines. The objec- 
tive of the cost study is to establish and demonstrate methods by which cost reductions 
may be achieved. Such reductions could significantly broaden the market for turbojet 
engines in general aviation. For these applications a control must be fully automatic and 
include all of the basic control functions. These include automatic engine startup, gov- 
erning with maximum speed limiting, acceleration limiting, and deceleration limiting. A 
discussion of these factors is presented in references 1 and 2. 

To satisfy the goal of low-cost, fuel- control system hardware should be less complex 
than that heretofore required for more conventional systems. Complex cams, trans- 
ducers, and servosystems should be avoided. The fuel-control concept which has evolved 
in the study can be mechanized to meet this criterion. The concept is based on the 
approximation of engine acceleration and deceleration limits with functions based on cor- 



rected parameters. The corrected parameter schedules are mechanized with hydraulic 
circuitry to form the control. 

At the time of the evaluation of the fuel- control concept, both an analog simulation 
of the J85-13 engine and an actual J85-13 engine were available for the program. The 
fuel- control evaluation was thus performed with the J85-13 simulation and engine. Al- 
though this engine is equipped with an afterburner, it is otherwise similar to the engine 
types considered in low-cost applications. During the test program, the afterburner was 
not used, and the compressor variable geometry and the variable exhaust nozzle were 
programmed with their normal operating schedules. 

The initial evaluation study was performed on an analog computer with a simulation 
of the J85-13 engine. The evaluation was then performed with a simulated fuel- control 
in combination with an electrohydraulic flow control and an actual J85-13 engine. 


CONTROL PRINCIPLES 


The fuel- control concept is based on the assumption that the maximum and minimum 
fuel limits can be described by linear functions of corrected engine parameters. A cor- 
rected parameter map illustrating these functions is presented in figure 1. The acceler- 
ation limit is defined by the relation 




( 1 ) 



Figure 1. - Corrected parameter map employed in 
experimental fuel control. 
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and the deceleration limit by the relation 




(2) 


(Symbols are defined in appendix A. ) These data were derived from approximations of 
the engine manufacturer's characteristics. As is shown in figure 2, different accelera- 



Figure 2. - Variation of J85-13 acceleration limits at different com- 
pressor inlet temperatures. 


tion limits are required as the inlet temperature varies to prevent excessive turbine inlet 
temperatures. As can be seen from equation (1) the effects of inlet temperature are 
neglected in the control. The lower acceleration limit is thus used for all inlet tempera- 
tures. With this simplification, the maximum allowable fuel cannot be used for the lower 
temperature range. It is expected that the resulting lower acceleration capability is ac- 
ceptable for general aircraft. 


CONTROL MECHANIZATION 

The hydraulic circuit employed in the fuel control is illustrated in figure 3. The cir- 
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Figure 3. - Schematic drawing of basic hydraulic circuit used in experimental fuel control. 


cuit generates an output flow rate W,p in accordance with the relation 

W t =M(A 1+ A v ) (3) 

A 1 

The pressure regulator bypasses excess flow from the main fuel pump to regulate the 
pump discharge pressure P to the value of P' . Thus, the pressure difference across 
the control pump is zero, and the pressure differences across orifices Aj and A v are 
equal. These pressure relations yield the following expressions for fuel flow: 

Wj_ - dN 


and 


W = 


A 

_v 

Ai 


W, 


The total fuel flow is 

W T - Wj + W y 

Substituting the expressions for and W y yields equation (3). 

It can be observed that the output flow (eq. (3)) is independent of discharge pressure 
P^ whenever the pressure regulator maintains a zero pressure drop across the control 
pump. In order to match the parameter map (see fig. 1) with the actual circuit (fig. 3), 
the bypass orifice Ay is composed of a series- parallel combination of orifices as shown 
in figure 4. 
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Control pump 




A mechanical schematic of the fuel control variable orifice system is shown in fig- 
ure 5. The areas of the various metering orifices are fixed by the spool positions which 
are, in turn, determined by the balance between the spring and pressure forces across 
the spool. The area of the governor orifices Ag thus depends on the net force across 
the first spool. The forces acting on this piston consist of the differential pressure 
(P - P, ) and the spring forces. The net forces applied are such that the governor area 
is closed at zero throttle position and a mechanical stop limits the maximum spool travel. 
The orifice areas A^ and Ag are fixed by the position of the second spool, which, in 
turn, is determined by the force balance between the compressor inlet pressure Pg and 


L 







the spool return springs. The orifice areas Ag and Ag are similarly determined by a 
force balance between the pressure rise across the compressor and the spool return 
springs . 

When the throttle is advanced, the force to the governor spool through the throttle 
spring is increased, thus advancing the governor spool until the force is balanced by the 
return spring and the differential pressure across the spool. For large throttle advances, 
the spool travel is limited by a mechanical stop. As engine speed increases, compressor 
discharge pressure Pg also increases to advance the third spool. This Spool opens ori- 
fice areas Ag and Ag to supply more fuel to the combustor. Concurrently, the differ- 
ential pressure (P a - P b> rises with speed thus applying a larger differential pressure 
across the orifice system and increasing the flow. The increase in (P - P^) also results 
in an increase in the return force on the governor spool and a corresponding reduction in 
Ag. An equilibrium condition is obtained whenever the power lever force applied to the 
springs is balanced by the force developed by the differential pressure (P - P b ) at a con- 
stant or equilibrium speed. 

ANALYSIS 

Equations relating the fuel- control flow rate to throttle position and engine param- 
eters such as speed, inlet pressure, and compressor pressure rise are developed in ap- 
pendix B. It is shown in appendix B that the total flow to the engine is given by 



The orifice area expressions are developed in appendix B as 

A C1 = A 3 + A 4 = c i + C 2 P 2 + C 3^ P 3 " *V 

A C2 " A 5 + A 6 = C 4 + C 5 P 2 + C 6^ P 3 " 

The governor orifice area Ag is a function of the spool position X g 

A 2 =f(X s ) (7) 

and the spool position X is a function of the throttle angle and engine speed 

b 
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( 8 ) 


g(a,N 2 ) - CrjOt - C 8 N 2 - C 9 
X s =g(a,N 2 ) g^X g 

X s= X s S >X s 


(9) 


Equations (4) to (9) define the behavior of the fuel control, relating the fuel-flow rate to 
engine speed, throttle position, and compressor inlet and discharge pressures. Because 
of the algebraic complexity of the set and the interaction of the control with the engine, 
the details of system operation will be discussed in conjunction with the system simulation. 

The engine acceleration and deceleration limits, however, can be visualized directly 
from equation (4). During rapid accelerations, the governor orifice area Ag becomes 
Ag, where A 2 is the full open governor orifice area. Therefore equation (4) becomes 


dN 

T (max) . 

A 1 



( 10 ) 


Because A^ and A^ are both functions of compressor inlet pressure and pressure 
rise, the acceleration limit is a function of speed, compressor inlet pressure, and com- 
pressor pressure rise. 

During deceleration, the governor orifice closes and equation (4) becomes 


W 


T(min) 


dN 


(A 


1 + A C2 ) 


(ID 


The deceleration limit is thus a function of speed, compressor inlet pressure, and com- 
pressor pressure rise, as defined. 

The analysis presented in appendix C relates the fuel- control design to engine fuel 
requirements. It is shown that the orifice areas A^ and A^ can be determined by 
considering expression (4) and the minimum and maximum engine fuel-flow requirements 
given by the linear relations 
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I 



(minimum flow) 


(12a) 


W n 


.m. = c,3 


N 

Vo 


'10 


w n 


-2VJL = c„ — + c 


N 

Vo 


'11 


12 


(maximum flow) 


(12b) 


Equations (12) are rewritten in appendix C to yield expressions for uncorrected minimum 
and maximum fuel flow: 


«T(min) = f C 10 P 3 
r 


W n 


— [( 
TD L 


T(max) L C 10 P 3 + C 12 P 2 + ( C ll " '“'10 


>> p 3 ] 


(13) 


— 2 2 

Comparing equations (13) with equations (10) and (11) and assuming that A 2 » re- 

sults in 


A d (Fi2 p 2 + (c n " c io) p a] 

l (14 

A C2 C 10 P 3 " A 1 

The equations for A^ and A^ 2 are expressed in terms of compressor inlet and dis- 
charge pressures. The coefficients C 1Q , C^, and C 12 can be determined directly 
from the engine fuel limit requirements. The acceleration schedule is thus based on ap- 
proximating the normal engine schedule to one based on equation (12b). 


PROCEDURE 

System Simulation 

A nonlinear, high-frequency simulation of the J85-13 engine (ref. 3) was used. In 
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this simulation, each compressor stage is represented by its pressure and temperature 
maps. The momentum equation is used to establish the flow dynamics from the stage 
pressures. Energy and mass balance equations and an equation of state are used to com- 
pute the pressures and temperatures of each stage volume. Engine speed and compres- 
sor inlet pressure and temperature are input variables, thus allowing for variable oper- 
ating conditions. An important feature of this dynamic model is that the compressor 
stall line is generated. 

The basic energy, continuity, and momentum equations are used for the combustor 
and turbine simulation. The turbine pressure and enthalpy relations are obtained with 
two-dimensional maps. An enthalpy balance between the compressor and turbine is com- 
puted to obtain shaft torque, from which engine acceleration and speed are generated. 

The normal variable geometry schedules (inlet guide vanes, interstage bleed, and ex- 
haust nozzle area) are included in the simulation. 


Simulated Fuel Control and Engine 

The fuel control was simulated on the analog computer and operated in conjunction 
with the J85-13 engine simulation. The information flow diagram of the control is pre- 



Figure 6. - Block diagram of fuel control. 
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sented in figure 6. The control simulation uses the basic component equations developed 
in appendix B with some modifications to avoid algebraic loops. 

Standard orifice flow equations are used to compute the flow across the system ori- 
fices. The pressure in chamber C (see fig. 4) is computed from the relation 

=-£- (W« - W„ - W 4 ) (15) 

c pv c 2 3 4 

The pressure in the volume between the fuel- control burner spray nozzles is computed 
from the expression 


< 16 > 

P V b 

The flow areas A g , A 4 , A g , and A g were simulated using equations (B18) to (B21). 

The equation for the governor spool position as a function of the spool forces is developed 
in appendix B as 


< p b - p a> A p - F b + - x s> - k 2 x s = + Dx s < 17 > 

Simulated Fuel Control with an Actual Engine 

The fuel control was simulated on a small analog computer and used in combination 
with an electro-hydraulic flow regulator to control an actual J85-13 engine. The engine 
was equipped with the original fuel control, exhaust nozzle area, variable geometry, and 
interstage bleed controls. Provisions are included, however, to transfer control action 
to research fuel controls and to operate the controls off schedule. 

The simulated fuel-control output signal was connected to a high-response electro- 
hydraulic flow regulator (ref. 4). Speed and pressure signals were fed back to the con- 
trol from engine mounted transducers. 


RESULTS AND DISCUSSION 

Performance studies of the fuel control were conducted both open loop and with en- 
gine feedback. The objective of this effort was to establish the steady- state and dynamic 
performance of the fuel control. 

The maximum and minimum values of fuel flow to speed parameter W t /NS as a 
function of pressure ratio supplied by the control are presented in figure 7. These fuel- 
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Figure 7 . - Acceleration and deceleration limits for experimental fuel control. 


flow limits correspond to equations (C17) and (C18) and approximately to equations (12a) 
and 12(b). As discussed in appendix C, the control mechanization is based on a design 
representing equations (C17) and (C18) and therefore includes a term inversely propor- 
tional to inlet pressure Pg. This term accounts for the spread observed for the mini- 
mum fuel flow limit in figure 7. 

The J85-13 simulation was used to demonstrate the control concept and to observe the 
steady-state and dynamic performance. The parametric curves of W T /N5 as a function 
of pressure ratio and throttle position are presented for sea-level conditions in figure 8(a). 
Acceleration occurs along a line of constant throttle position (X ff line) to the final speed. 
It can be observed that at the extremes, the X ff lines converge to maximum and mini- 
mum limits. Similar curves for altitude conditions (25 000 ft (7600 m), M = 0.65) are 
presented in figure 8(b). It can be observed that the minimum Wp /N6 limit varies as 
a function of inlet pressure P2, as given by equation (Cl 7). In addition, the Xq, lines 
generated by the control are not identical to the sea-level condition. 

Curves of fuel flow to pressure parameter Wrp/Pg as a function of corrected engine 
speed and throttle position are shown for sea-level conditions in figure 9(a). With these 
coordinates, the slopes of the constant lines between the maximum and minimum 
fuel limits are related to the gain of the proportional governor. As indicated by figure 
9(a), in the vicinity of the acceleration limit, these lines are nearly horizontal because 
the governor is overridden by the acceleration limit. As the operating condition is ap- 
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Figure 10. - Comparison of acceleration schedules 


proached, governor action is restored, and the gain increases to provide good control 
performance. Similar curves for altitude condition are presented in figure 9(b). 

Satisfactory fuel- control performance requires that the maximum fuel limit restrict 
engine operation to safe operating regions. Figure 10 presents the acceleration schedules 
for the experimental and engine manufacturer’s controls. The engine steady-state oper- 
ating lines are included to indicate the available acceleration margin. In addition, typical 
acceleration histories from near 70 to 100 percent corrected speed are shown for the ex- 
perimental control. As the corrected speed begins to approach the 100-percent value, 
the trajectory leaves the acceleration schedule, and normal proportional control is re- 
stored. It can be seen from the tabulations of representative response times shown in 
table I, that for sea-level conditions the experimental and normal engine controls have 
equivalent response times. The experimental control regulates within the boundaries of 
the engine manufacturer's control except for a small excursion near idle speed. 

The variation of mechanical engine speed as a function of equivalent throttle displace- 
ment X in steady state is presented in figure 11 for sea-level and altitude conditions. 

If the altitude dependency noted here is undesirable, inlet pressure or temperature cor- 
rections may be included in the design. The indicated values give the actual linear dis- 
placement of the governor spring and must be converted through a cam. It can be con- 
cluded from these curves that a physical limit must be designed into the power lever so 
that the engine design speed will not be exceeded. As an alternate approach, an altitude 
compensated throttle can be designed for the control. 
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TABLE I. - SUMMARY OF ENGINE RESPONSE TIME 


[Sea level static. ] 


Speed , 
percent 

Response time, sec 

J85-13 engine 
with normal 
control 

J85-13 engine 
with simulated 
experimental 
control 

Simulated 
engine with 
simulated 
experimental 
control 

70 - 100 

0.71 

0. 7 

0.80 

80 - 90 

. 5 

. 3 

. 34 

80 - 95 

. 5 

.45 

.47 

80 - 100 

a . 5 

.5 

.48 

90 - 100 

a . 3 

— 

.30 


initial time to desired value additional settling time 
about 1. 5 sec. 



40 ' 

0 


.10 . 20 . 30 . 40 . 50 . 60 

Equivalent throttle displacement, in. 


0 .2 .4 .6 .8 1.0 1.2 1.4 

Equivalent throttle displacement, cm 

Figure 1L - Equivalent throttle displacement as function of actual 
speed. 




(a) Acceleration, 70 to 100 percent speed. 


Figure 12. - Dynamic response of simulated engine to throttle command (sea level) for 
engine acceleration of 70 to 100 percent 

The J85-13 engine and fuel-control simulations were also subjected to throttle dis- 
turbances to determine the acceleration and deceleration characteristics at sea level and 
at 40 000 feet (12 200 m) altitude (M = 0.91). Figure 12(a) presents a recording of the 
engine acceleration from 70 to 100 percent speed. For this throttle advance, the gover- 
nor area was forced fully open and the engine accelerated to the desired operating speed 
in about 0. 8 second. Engine acceleration response from 90 to 100 percent is illustrated 
in figure 12(b). For this case, the magnitude of the disturbance was insufficient to force 
the governor area to the fully open position. The response time for the 90 to 100 percent 
speed change was 0. 3 second. 

Engine deceleration from 100 to 70 percent is illustrated in figure 13. This throttle 
chop caused the governor orifice to close. The time required to decrease the speed to 
the desired condition is about 1.6 seconds. 

The system performance was also evaluated at 40 000 feet (12 200 m) and Mach 0.91. 
A representative recording illustrating acceleration from 90 to 100 percent speed is pre- 
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Compressor Fuel flow, Speed, percent Throttle position 

discharge Ib/hr (kg/hr) 

pressure, 




Compressor Fuel flow, Speed, percent Throttle position 
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pressure, 




Figure 13. - Dynamic response of simulated engine to throttle chop at sea level. Engine 
deceleration, 100 to 70 percent speed. 



sented in figure 14. The time required to achieve the steady-state condition is about 
3.0 seconds. A similar response for engine deceleration is shown in figure 15. The re- 
sults indicate that a 100 to 90 percent speed change was accomplished within 3. 0 seconds. 

The response of the actual J85-13 engine with the simulated control is shown in 
figure 16. These tests demonstrated the effectiveness of the acceleration limit in pre- 
venting engine stall at sea-level conditions. A typical response from 70 to 99 percent 
speed is illustrated in figure 16(a). This speed change was accomplished within 
0. 7 second. A similar recording showing a 90 to 99 percent speed change is presented 
in figure 16(b). This step change was achieved within 0. 3 second. During the test pro- 
gram with the actual J85-13 engine, step commands were restricted to 99 percent 


a> 




a> 







28.76 (19.8) 


0. 2 sec ■ 

I II I 


Figure 14. - Dynamic response simulated engine to throttle command at 40 000 feet 
(12 200 m) and Mach 0.91. Engine acceleration, 90 to 100 percent speed. 
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(b) Acceleration, 90 to 99 percent speed. 
Figure 16. - Concluded. 


maximum speed or less. The restriction was imposed to eliminate any possibility of af- 
terburner light-off. 

During engine acceleration, the exhaust nozzle was incorrectly scheduled. This 
schedule resulted in lower steady-state fuel flow. For this reason, the results obtained 
from the engine simulation at sea level (fig. 12) do not fully correspond to the actual 
engine tests (fig. 16). The exhaust nozzle area variation of the simulated engine was 
matched to that of the test engine. The performance for a 70 to 99 percent speed dis- 
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Figure 17. - Dynamic response of simulated engine with off-schedule 
exhaust nozzle area at sea level. Engine acceleration, 70 to 99 per- 
cent speed. 


turbance is shown in figure 17. A settling time of 0. 8 second was obtained. The 
2-percent speed overshoot can be attributed to the difference in exhaust nozzle movement 
A similar recording for a 90 to 99 percent speed change is presented in figure 18. This 
step change was accomplished within 0.4 second. A tabulation of representative system 
response times is presented in table I. The tabulation includes response times obtained 
from the actual engine and its normal controls, the actual engine with the simulated ex- 
perimental control, and the engine simulation with the simulated control. The data pre- 
sented indicate that the performance obtained from the control concept examined in this 
report compares favorably to that obtained with the original engine controls. 


23 





Figure 18. - Dynamic response of simulated engine with off-schedule ex- 
haust nozzle area at sea level. Engine acceleration, 90 to 99 percent 
speed. 


CONCLUSIONS 

The fuel- control tests with the actual J85-13 engine indicate that satisfactory steady- 
state and dynamic performance can be achieved. The acceleration limiter was effective 
in preventing engine stall. 

The tests conducted with the simulated J85-13 engine at sea-level and altitude condi- *• 

tions indicate satisfactory performance over the desired operating range. The control 
effectively restricted engine operation within permissible limits. 
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The results of the analytical and test programs indicate that the techniques employed 
in the design of the generalized parameter hydraulic fuel control are satisfactory for 
single-spool turbojet engines. The original objectives of the control design were met, 
and the results show that the inherent simplifications do not affect the dynamic perform- 
ance of the control- engine system. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 4, 1969, 

720-03. 
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APPENDIX A 


SYMBOLS 

A C1 sum of orifice areas A g and A 4 , in. 2 ; cm 2 
A C2 sum of orifice areas A 5 and A g , in. 2 ; cm 2 

O p 

A^ general orifice area, in. ; cm 

p p 

A n spray bar nozzle area, in. ; cm 

p p 

Ap governor piston area, in. ; cm 

p p 

Ay bypass orifice area, in. ; cm 

p p 

Aj reference orifice area, in. ; cm 

p p 

A 2 governor area, in. ; cm 

A 2 maximum governor area, in. 2 ; cm 2 

A 3 orifice area, f(P 3 , P 2 ), in. 2 ; cm 2 

A 4 orifice area, f(P 2 ), in. 2 ; cm 2 

Ag orifice area, f(P g , P 2 ), in. 2 ; cm 2 

Ag orifice area, f(P 2 ), in. 2 ; cm 2 

a 34 area coefficient 

a 32 area coefficient 

a 4 j area coefficient 

a 42 area coefficient 

a g4 area coefficient 

a g2 area coefficient 

a g j area coefficient 

a g2 area coefficient 

C i constants (i = 1, . . . , n) 

D damping coefficient, (lbf)(sec)/in. ; (N)(sec)/m 
d pump flow coefficient, (lb) (hr)/% speed; (kg) (hr)/% speed 
F b bias force, lbf; N 
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k 


2 

M 

N 

P a 

p ; 

p b 

p b 


c 

A Pi 


s 

T 2 

t 

V b 

V c 

AWj 

AW, 


W i 


W 


N 


W n 


w 


w 


T (max) 
T (min) 


W. 


p p 

orifice flow conversion constant, lbm/(hr)(in. ) (lbf/in. ) 
(N/cm 2 ) 1 ' 2 

spring constant, lbf/in.; N/cm 

spring constant, lbf/in.; N/cm 

mass of governor spool, lbm; kg 

engine speed, percent of rated speed 
2 2 

pressure, lbf/in. ; N/cm 

p p 

pressure, lbf/in. ; N/cm 

p p 

pressure, lbf/in. ; N/cm 

2 2 

rate of change of pressure, (lbf)(in. )/sec; N/cm /sec 
2 2 

pressure, lbf/in. ; N/cm 

2 o r 

rate of change of pressure, (lbf/in. )/sec; (N/cm )/'sec 

2 2 

general differential pressure, lbf/in. ; N/cm 

p p 

reference pressure, lbf/in. ; N/em^ 

2 2 

compressor inlet pressure, lbf/in. ; N/cm 

2 2 

compressor discharge pressure, lbf/in. ; N/cm 

p p 

burner pressure, lbf/in. ; N/cm 
Laplace operator 
inlet temperature, °R; K 
time, sec 

3 3 

volume of chamber b, in. ; cm 

3 3 

volume of chamber c, in. ; cm 
net fuel flow, lbm /hr; kg/hr 
net fuel flow, lbm/hr; kg/hr 
general fuel flow, lbm/hr; kg/hr 
spray bar fuel flow, lbm/hr; kg/hr 
total fuel flow, lbm/hr; kg/hr 
acceleration fuel flow, lbm/hr; kg/hr 
deceleration fuel flow, lbm/hr; kg/hr 
bypass fuel flow, lbm/hr; kg/hr 


1/2 


kg/ (hr) (cm 2 ) 


27 


fuel flow through orifice Ap lbm/hr; kg/hr 
Wg fuel flow through orifice A 2 , lbm/hr; kg/hr 

Wg fuel flow through orifice Ag, lbm/hr; kg/hr 

fuel flow through orifice A^, lbm/hr; kg/hr 
Wg fuel flow through orifice Ag, lbm/hr; kg/hr 

Wg fuel flow through orifice A g , lbm/hr; kg/hr 

X ff equivalent linear displacement of throttle, in. ; cm 
X g linear displacement of governor spool, in. ; cm 
X 0 maximum displacement of governor spool, in. : cm 

fa 

a angular throttle motion, rad 

Ace change of throttle motion, rad 

(3 bulk modulus, lbf/in.^; N/cm^ 

6 P 2 /14.7; P 2 /10. 1 

yfe JT 2 /518.1; ^/T 2 /288 

o o 

p fuel density, lbm/in. ; kg/cm 
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APPENDIX B 


SYSTEM EQUATIONS 

The weight flow across the various orifices of the fuel control can be described by 
equations of the form 

Wj = kAj yXpT (Bl) 

where the density term is assumed constant and included in the coefficient k. For the 
six orifices of the system shown in the schematic of figure 1, this equation becomes 


3 

h-* 

II 

£ 

h-* 

»/ P a- 


(B2) 

W 2 =kA 2 , 

^a- 

P c 

(B3) 

3 

CO 

II 

% 

CO 

^c- 

p b 

(B4) 

W 4 &kA 4’ 

V^c- 

p b 

(B5) 

W 5 - kA 5 

V P a- 

p b 

(B6) 

W 6 =kA 6- 

V P a- 

p b 

(B7) 


and the final fuel flow through the spray bar nozzle into the engine is 

W N = ka n (B8) 

Since the flow through the governor orifice A 2 passes through Ag and A^, 

W 2 = W 3 + W 4 (B9) 

Equations (B3) to (B5) can be combined with equation (B9) to give an expression for the 
intermediate pressure P c : 

p A|P a - (A 3 + A 4 ) 2 P b 

c 9 9 

(Ag + A 4 ) + Ag 
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Defining as A^ 


^ p a + A Sl p b 

A^ + A 2 
A 2 + A C1 


(BIO) 


and 


w 3 + w 4 = k[ 


A 2 A 2 
A C1 A 2 

kA Cl + A 2/ 


1/2 


( P a - P b> 


1/2 


(Bll) 


If the flows Wg and Wg are summed 


W 5 + W 6 =kA C2( P a- P b> 


1/2 


(B12) 


where 


A C2 “ A 5 + A 6 

The flow through the control pump is proportional to the engine speed, hence 


Wj = dN 


(B13) 


Combining equations (B2) and (B13) yields 

a p b 


(P Q - P h ) 1/2 = — 

kA, 


(B14) 


Substituting the differential pressure from equation (B14) into equations (Bll) and (B12) 
results in 


W3 + W4 - 


2 2 l/ 2 

A C1 A 2 dN 


A 2 + A 2 A 1 

A C1 A 2 1 


(B15) 


+ w 6 = a c ® 
A 1 


(B16) 


30 



The total fuel flow to the engine can now be computed from equations (B13), (B15), and 
(B16) as 


A i 



(B17) 


Equation (B17) relates the engine fuel flow, shaft speed, and the various area terms. 
The flow areas are determined as illustrated in figure 4, by the force balance between 
the pressure and spring forces applied across the spools of the fuel control. 

The areas and Ag, which are determined by the position of the second spool, 
can thus be described by 


A 4 - a 41 + a 42 P 2 

(B18) 

A 6 = a 61 + a 62 P 2 

(B19) 


Similarly, Ag and A^ are determined by the position of the third spool: 


A 3 “ a 31 + a 32^ P 3 " P 2^ 
Ag = a 51 + a 52 (P 3 - P 2 ) 


(B20) 

(B21) 


The combined area terms A^ and A^g of equation (B17) can be formed from equa- 
tions (B18) to (B21), where 


A C1 _ A 3 + A 4 ~ ^1 + C 2 P 2 + C 3^ P 3 P 2^ (B22) 

A C2 = A 5 + A 6 = C 4 + C 5 P 2 + C 6^ P 3 " P 2^ (B23) 

The coefficients C- can be determined from the areas, spring constants, and flow area 
against stroke characteristic of the system. The area characteristics for the orifices 
controlled by the second and third spools are linear functions of spool position, but the 
governor orifice area is a nonlinear function of the stroke of its spool. 

It can be seen from figure 5 that, as the throttle is advanced through an angle Aa, 
the throttle cam converts the rotation to a linear advance against the throttle spring. 

The spring, in turn, raises the forces on the spool tending to increase the governor ori- 
fice area Ag. If the throttle advance Acn is sufficiently large, the spool motion is 

limited by a mechanical stop. The spool response to throttle position can be described 
by 
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< P b - P a> A p - F b + k l «a - x s> - Vs = ^s + Dx s < B24 > 

Neglecting the dynamics of the mass spring system reduces equation (B24) to 

< P b - P a> A p - F b + k l ^ - x s> - Vs = 0 < B25 > 

Since 

(P a - P b ) « N 2 


The equation reduces to 


-N 2 A, 


+ k l X a 


- 0^ + k 2 )X g - F b = 0 


(B26) 


A nonlinear relation could exist between power lever and throttle spring position. How- 
ever, for this analysis, this function was assumed to be linear. Equation (B26) becomes 


g(a,N 2 ) - C 7 a - C g N 2 - C 9 
X s =g(a,N 2 ) gs=Xg 


X 


s 



> X„ 


(B27) 


The governor orifice area A 2 is a logarithmic function of the throttle position X g : 


A 2 =f(X g ) (B28) 

The nonlinear characteristic was selected to compensate for the higher ratio of engine 
speed to fuel gain encountered for the lower speed region. By this method the overall 
system gain can be modified to prevent instability. 
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APPENDIX C 


FUEL CONTROL ORIFICE SIZING 


Once the basic design of the fuel control is established, the control components must 
be sized to satisfy the fuel requirements of the engine. An expression for the total flow 
from the control to the engine was developed in appendix B as 


W 


T " 



A 1 + A C2 


| A C1 A 2 

V A C1 + A 2, 


(Cl) 


For minimum flow the governor area A 2 is fully closed and expression (Cl) re- 
duces to 


W T(min) 


(A 1 + A C2^ 

A 1 


(C2) 


Similarly, for maximum flow the governor area A 2 is fully open and equation (Cl) 
becomes 


w - Nd 

W T (max) 

A 1 


A 1 + A C2 


A C1 A 2 



where A 2 = A 2max 

n o 

Assuming that A| » A C1 , equation (C3) reduces to 


w T(max) a f 5 <A 1 + A ' 


C2 + A C1^ 


(C3) 


(C4) 


The fuel- control concept is based on the assumption that the minimum and maximum 
fuel-flow limits to corrected speed ratio can be expressed as a linear function of com- 
pressor pressure ratio. Specifically, it is assumed that 

w T 

— 7 ~ p 

6 v . g = Q, n — (minimum flow) (C5) 

JL 10 P 2 
yfe 


33 



(maximum flow) 


(C6) 


w T 

sVe 

N 

V? 


= c 


11 


+ c 


12 


Rewriting equations (C5) and (C6) yields 

W 


T(min) ~~~ ( C 10 P 3^ 


(C7) 


^T(max) [ C 10 P 3 + C 12 P 2 + (C 11 " C 10 )P s] (C8) 

r 

A comparison between the fuel-control equations (C2) and (C4) and the engine fuel 
requirements (C7) and (C8) yields the following expressions from which the orifice areas 
and A C2 may be computed 


[ C 12 P 2 + < C 11 - C 10> P 3 ] 

(C9) 

A 1 


A C2 = p H C 10 P 3 - A 1 
r 

(CIO) 


The coefficients C-^q, C-q, and C 12 must be evaluated from the minimum and maxi- 
mum engine fuel requirements. For the J85-13 engine considered in the analysis of this 
report, these requirements were approximated by selecting 


Cio = 1-54 (lbm/hr)/% speed 
= 0.70 (kg/hr )/% speed 


U. S. customary units 
SI units 


= 4. 15 (lbm/hr)/% speed 
= 1.88 (kg/hr )/% speed 


U. S. customary units 
SI units 


Cj2 = 8. 195 (lbm/hr)/% speed 
= 3. 72 (kg/hr )/% speed 


U. S. customary units 
SI units 


Substituting these values into the area relations of equations (C9) and (CIO) and using 
values for d, P , V, and A^ from table II results in 
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TABLE H. - CONSTANTS 


2 2 

Reference orifice area, A^, in. ; cm 

3.24X10" 2 ; 20.9X10 -2 

— 2 2 

Maximum governor area, A2, in. ; cm 

0.15; 0.968 

2 2 

Governor piston area, A^, in. ; cm 

1; 6.45 

Conversion constants: 


Cp in. 2 ; cm 2 

2.2X10' 2 ; 14.2X10 -2 

C 2 , in. 2 /(lbf/in. 2 ); cm 2 /(N/cm 2 ) 

1.47X10' 2 ; 13.8X10" 2 

C,, in. 2 /(lbf/in. 2 ); cm 2 /(N/cm 2 ) 

3.55X10" 4 ; 33.2X10' 4 

C4, in. ; cm* 5 

1.22X10" 2 ; 7.87X10' 2 

C 5 , in. 2 /(lbf/in. 2 ); cm 2 /(N/cm 2 ) 

2. 10X10 -4 ; 19.7X10' 4 

Cg, in. 2 /(lbf/in. 2 ); cm 2 /(N/cm 2 ) 

2.10X10" 4 ; 19.7X10' 4 

C 10 , (lbm/hr)/% speed; (kg/hr )/% speed 

1.54; 0.70 

C n , (lbm/hr)/% speed; (kg/hr )/% speed 

4. 15; 1. 88 

C 12 i (lbm/hr)/% speed; (kg/hr)/% speed 

8.195; 3.72 

Damping coefficient, D, (lbf)(sec)/in. ; N-sec/cm 

0. 18; 0. 315 

Pump flow coefficient, d, (lbm/hr)/% speed; (kg/hr)/% speed 

1.62; 0.737 

Bias force, F^, lbf; N 

5; 22.24 

Orifice flow conversion constant, k, (lbm/hr)in. 2 (lbf/in. 2 ) 4//2 ; 

1.02X10 4 ; 863 

(kg/hr) cm^N/cm 2 ) 1 ^ 2 


Spring constant, kp lbf/in. ; N/cm 

75; 131.3 

Spring constant, k 2 , lbf/in. ; N/cm 

25; 43.8 

Mass of governor spool, M, lbm; kg 

0. 327xl0“ 2 ; 0.0574 

Reference pressure, P , lbf/in. 2 ; N/cm 2 

14.7; 10.1 

Maximum linear displacement of governor spool, X g , in. ; cm 

0.20; 0.508 

Volume of chamber b, V, , in. 2 ; cm 2 

6; 98.32 

Volume of chamber c, V c , in. ; cm 

20; 327.7 

Bulk modulus, / 3 , lbf/in. 2 ; N/cm 2 

1.5X10 5 ; 1 . 034x1 0 2 

Fuel density, p, lbm/in. 2 ; kg/cm 2 

0.0291; 0. 806X10" 2 


A C1 = 3. 551X10" 4 (P 3 - P 2 ) + 14. 7xl0" 4 P 2 
= 33.2X10" 4 (P 3 - P 2 ) + 138x10 4 P 2 


U. S. customary imits 
SI units 


i 


(Cll) 


•,-4 


•^C2 ~ 2.095x10 ** (Pg - P 2 ) + 2.095x10 4 P 2 - 3.24x10 3 U. S. customary units 
= 19. 7X10 " 4 (P 3 - P 2 ) + 19. 7x10 " 4 P 2 - 20. 9X10' 3 SI units 


-3 


(C12) 


Equations (Cll) and (C12) show that the equivalent areas A C1 and A C2 can be formed 
from the combination of an area proportional to compressor pressure rise and an area 
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proportional to compressor inlet pressure. These expressions can be used with equa- 
tions (B23) and (B24) of appendix B to compute the following fuel control orifice areas: 

Aqj = Ag + (Cl 3) 

where 

Ag = 3. 55x10 (Pg - Pg) + a-gj U.S. customary units 

= 33. 2x10 -4 (Pg - Pg) + ag^ SI units 

-4 

A^ = 14. 7x10 Pg + a^ U.S. customary units 

= 138x10 ~ 4 ^2 + a 41 SI un ^ s 

^C2 = ^5 + -^6 (C14) 

where 

Ag = 2. 095x10 _4 Pg - 3.24X10 -3 U.S. customary units 
= 19. 7X10 -4 P 2 - 20. 9X10" 3 SI units 

-4 

Ag = 2.095x10 (Pg - Pg) + a^ U.S. customary units 

= 19. 7xl0~ 4 (Pg - Pg) + a^ SI units 

In the simulation of the actual fuel control, the area expressions (Cll) and (C12) were 
not exactly duplicated. The corresponding expressions for the design are: 

A C1 = 3. 55X10 -4 (Pg - Pg) + 14. 7xl0 -4 Pg - 22X10" 4 
= 33.2X10" 4 (P 3 - Pg) + 138x10 “ 4 Pg - 141. 9x10 " 4 

A c2 = 2. 10X10' 4 (Pg - P 2 ) + 2. 10X10 -4 Pg - 12. 2xl0 -4 
= 19. 7X10' 4 (Pg - Pg) + 19. 7x10 “ 4 Pg - 78. 7xl0 -4 
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Equations (C15) and (C16) are equivalent to the following engine minimum and maximum 
fuel-flow requirements: 


For minimum flow 


W n 


M. = 1.544 12. + iLM 
N 

V® 


P 2 P 2 


= 0.7l3 + i^3 
P 2 P 2 




U.S. customary units 


SI units 


(C17) 


For maximum flow 


Wr 


5y fe _ 

N 


4.153—3 + 8.195 - 1 - 2 - 


1.88—3+ 3.72 - 


U. S. customary units 


> (C 18 ) 


SI units 


It can be seen that equations (C15) and (C16) are close to equations (Cll) and (C12). 
Equations (C17) and (C18) are close to equations (C5) and (C6) except that a term in- 
versely proportional to inlet pressure Pg is included. The effect of this term is illus- 
trated by the scatter shown in figure 7. Since the purpose of the simulation was to dem- 
onstrate the feasibility of the design, the analytical work of this report is based on equa- 
tions (Cl 5) and (C16). 
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